DIN SRC. 4 Ss} losy 


Jun \oe| 


DAVID W. TAYLOR NAVAL SHIP 
RESEARCH AND DEVELOPMENT CENTER 


Bethesda, Maryland 20084 


DTNSRDC-81/054 


PERIODIC SINGLE-BLADE LOADS ON PROPELLERS 
IN TANGENTIAL AND LONGITUDINAL WAKES 


WHO] 


DOCUMENT uy 
COLLECTION 


Robert J. Boswell 
Stuart D. Jessup 
Ki-Han Kim 


APPROVED FOR PUBLIC RELEASE: DISTRIBUTION UNLIMITED 


Presented at the 
PROPELLERS ’81 SYMPOSIUM 
The Society of Naval Architects and Marine Engineers 
Virginia Beach, Virginia, May 26—27, 1981 


INGLE-BLADE LOADS ON PROPELLERS IN TANGENTIAL 


wn 
ea SHIP PERFORMANCE DEPARTMENT 
= RESEARCH AND DEVELOPMENT REPORT 
— 
< 
2 
2) 
=) 
= 
2 =e 
ce 
| June 1981 DTNSRDC-81/054 
D3 | | 


20 91/054 


—e—— 
lerinersedes 3960/44) 


MAJOR DTNSRDC ORGANIZATIONAL COMPONENTS 


DTNSRDC 
COMMANDER 00 
TECHNICAL DI lar 


OFFICER-IN-CHARGE 
CARDEROCK 


OFFICER-IN-CHARGE 
ANNAPOLIS 


SYSTEMS 
DEVELOPMENT 
DEPARTMENT 


AVIATION AND 
SURFACE EFFECTS 
DEPARTMENT 


SHIP PERFORMANCE 
DEPARTMENT 


15 


COMPUTATION, 
MATHEMATICS AND 
LOGISTICS. Per An eENy 


STRUCTURES 
DEPARTMENT 


SHIP ACOUSTICS PROPULSION AND 
DEPARTMENT AUXILIARY SYSTEMS ae 
DEPARTMENT 


] 
| 
| 


SHIP MATERIALS CENTRAL 


ENGINEERING INSTRUMENTATION 
DEPARTMENT Of DEPARTMENT 


I 


VIM 


0 0301 0032540 1 


MBL/WHOI 


MM 


NDW-DTNSRDC 3960/43 (Rev. 3-78) 
GPO 728 741 


UNCLASSIFIED 


LE 
SECUMITY CLASSIFICATION OF THIS PAGE (When Data Entered) 1 


READ INSTRUCTIONS 
REPORT DOCUMENTATION PAGE BEFORE COMPLETING FORM 
feiecoey Poe ee ae 


DINSRDC-81/054 
4. TITLE (and Subtitle) 5. TYPE OF REPORT & PERIOD COVERED 


PERIODIC SINGLE-BLADE LOADS ON PROPELLERS IN 
TANGENTIAL AND LONGITUDINAL WAKES 


Ship Performance Dept. 
Research and Development 


6. PERFORMING ORG. REPORT NUMBER 


8. CONTRACT OR GRANT NUMBER(s) 


. AUTHOR(s) 
Robert J. Boswell 


Stuart D. Jessup 
Ki-Han Kim 


10. PROGRAM ELE 
AREA & WORK 


NT, PROJECT, TASK 
IT NUMBERS 


. PERFORMING ORGANIZATION NAME AND ADDRESS 


David W. Taylor Naval Ship Research 
and Development Center 
Bethesda, Maryland 20084 
11. CONTROLLING OFFICE NAME AND ADDRESS 
Naval Sea Systems Command (05R) 

Ship Systems Research and Technology Group 
Washington, D.C. 20362 


14. MONITORING AGENCY NAME & ADDRESS(if different from Controlling Office) 


ME 
UN 


(See reverse side) 


12. REPORT DATE 
June 1981 
13. NUMBER OF PAGES 


88 


SECURITY CLASS. (of this report) 


15. 


Naval Sea Systems Command (524) 
Propulsion Line Shafting Equipment Division 
Washington, D.C. 20362 


UNCLASSIFIED 


1Sa. DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


. DISTRIBUTION STATEMENT (of this Repost) 


APPROVED FOR PUBLIC RELEASE: DISTRIBUTION UNLIMITED 


. DISTRIBUTION STATEMENT (of the abstract entered in Block 20, if different from Report) 


. SUPPLEMENTARY NOTES 


. KEY WORDS (Continue on reverse side if necessary and identify by block number) 
Marine Propeller Loads 
Propeller Research 


Propulsion 
Model Experiments 


. ABSTRACT (Continue on reverse side if necessary and identify by block number) 
Periodic single-blade loads were measured on model propellers in in- 
clined flow and in a single-cycle circumferentially non-uniform longitudinal 
velocity field. The experimental results were correlated with predictions 
by the following methods: 
1. A quasi-steady procedure developed by McCarthy at David W. Taylor 
Naval Ship Research and Development Center (DTNSRDC), 


(Continued on reverse side) 


DD ; FORM 1473 EDITION OF 1 NOV 65 1S OBSOLETE 


JAN 73 UNCLASSIFIED 


S/N 0102-LF-014-6601 eee 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


UNCLASSIFIED 


—__ eee 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


(Block 10) 


Task Area S0379-SLOO1 
Task 19977 
Work Units 1544-296 and 1544-350 


(Block 20 continued) 


2. An unsteady lifting surface theory developed by Tsakonas and 

his colleagues at Davidson Laboratory, 

3. An unsteady lifting surface theory developed by Kerwin and Lee 

at MIT, 

4, A refinement by Kerwin to the method of Kerwin and Lee to 

consider the inclination of the propeller slipstream. 

In inclined flow, all four of the calculation methods evaluated 
consistently underpredicted the experimental values of the periodic 
single-blade loads. The method of Kerwin, which considers the in- 
clination of the slipstream relative to the propeller axis, produced 
the best correlation with experimental values. These correlations 
show that the inclination of the propeller slipstream relative to the 
propeller axis significantly influence the periodic single-blade loads. 
The importance of this inclination increases with increasing time- 
average loading. 

In longitudinal flow, all of the calculation procedures predicted 
periodic single-blade axial force at design advance coefficient to 
within 20 percent of the experimental values, but agreement was not as 
good at off-design advance coefficients. The method of Kerwin and Lee 
produced the best overall correlation with experimental results 
considering both amplitude and phase of the periodic single-blade axial 
force over a range of advance coefficients. 


UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 


LIST OF FIGURES ... . 
LIST OF TABLES .... 
ABSTRACT 0 00010 60 
NOMENCLATURE 6G 010 
INTRODUCTION .... - 
BACKGROUND .... «© « 


EXPERIMENTAL TECHNIQUES 
WAKE SCREEN AND WAKE 


SURVEY 


TABLE OF 


EXPERIMENTAL CONDITIONS AND PROCEDURES 


,/DATA ACQUISITION AND ANALYSIS 


ACCURACY . . ... -» 


EXPERIMENTAL RESULTS . 
LOADING COMPONENTS . 


CENTRIFUGAL AND GRAVITATIONAL LOADS 


e 


CONTENTS 


VARIATION OF LOADS WITH ANGULAR POSITION. . . 
EFFECT OF PLATE CLEARANCE 


EFFECT OF ADVANCE COEFFICIENT 


CORRELATION BETWEEN EXPERIMENTAL RESULTS AND THEORETICAL PREDICTIONS 


THEORETICAL METHODS 


CORRELATIONS IN TANGENTIAL WAKES 
CORRELATIONS IN LONGITUDINAL WAKES 


SUMMARY AND CONCLUSIONS 


ACKNOWLEDGMENTS ... .~ 


REFERENCES ...-. -@ 


e 


iii 


Page 


10 


qa! 


12 


13 


LIST OF FIGURES 


Components of Blade Loading .......... 


Sign Convention for Ship Wake Velocity Components 


Distribution of Wakes in Propeller Disks ... 


Experimental Variation of Loading Components 
witheB laden Angular PO SistOmeny yeiiei iefetel memttelteatenite 


Effect of Plate Clearance on First Harmonic 
Loading in Tangential Wake ........2.-. © 


Experimental Variation of First Harmonic Blade 
Loads on Propeller 4661 with Advance Coefficient 


Comparison of First Harmonic Blade Loading 
Coefficients in Tangential Wakes ...... . 


Comparison of Various Experimental Blade 
Loading Coefficients in Longitudinal Wakes .. 


Amplitudes of First Harmonic Blade Loads in 
Tangential Wakes; Correlation Between Experiment 
E\igiel GNaVeoIRS? Gg 6) oO Go O65 6G 6 G6 00 06 6 8 


Phases of First Harmonic Blade Loads in Tangential 


Wakes; Correlation Between Experiment and Theory 


Variation of Periodic Blade Loads 
VAL Soeuee iheoulaligeneteiy od 6 6 6 6 6 6 0 6 0 0 


Amplitudes of First Harmonic Blade Loads in 
Longitudinal Wakes; Correlation Between 
EXperamenteandm THEORY / te) sve) 1s) Nell tel a tolacoMie Menor ol its 


Phases of First Harmonic Blade Loads in Longitudinal 


Wakes; Correlation Between Experiment and Theory 


iv 


Page 


12 


14 


18 


19 


19 


LIST OF TABLES 


1 - Characteristics of Propellers 4661, 4710, and 4402 .... 


2 - Comparison of Experiment and Theory in Inclined Flow - 
Amplitudes and Phases of First Harmonic Loading Coefficients 
AS DASICI Wo 6 o 6! 0) 0 6) oo O60 0 16 oO Go oN vol 6d) og ol 


3 - Comparison of Experiment and Theory in Longitudinal Flow - 
Amplitudes and Phases of First Harmonic Loading Coefficients 
aie Iyaeileiy Oo 666-6 1G oO OO Olkd WON Otomo G uOr dy oO) GO 1d 6 6) 6 


i 
Fal 
Hae 
1 " 
P i 
% rY ‘ i tos i 
Pan ee S oe 


f } rh r 
Lah : Tn \\ ER 

Pore 

aT ¢ ar ko oe 7 

HY % 
on ee it UNS ' 

es Paki 
Fie i “i My 

Sy, i H ’ f 
* 
hy rah ; 
“ x Ate 
uv : Dri 
a A 
} \ A Ls 
' 7 wh! W 
me) 
th \ Hos Oy 
" 
- i i} 


POW 


¢ od 
nee 


ii A Mahe sat Ai : 


ayy 


‘ m dies voa hg F ; 
f 4 v ri iV 4 4 
Me hie 04 Paw 
Lire Boe ry an eee ' 


BM ete \ Did F nines 


Fie a‘ cn ine ro 


ids nT 


ie ey ie ine vee e. 


: cata: "iss sa 1) aad 8 


a te wut ene duties 


sy 
n eye 


ar ile he steele as 
Ane Lae ane | Pay meneity, bas 


msi or a Oe vane 


¥ t 


Ee 


Ee |e 


5h RE oe 


Rew oe: Manet: 7 4 i “= 


a 


nvAM, Lethe: day” 


ABSTRACT 


Periodic single-blade loads were 
Measured on model propellers in inclined 
flow and in a single-cycle circumferen- 
tially non-uniform longitudinal velocity 
field. The experimental results were 
correlated with predictions by the 
following methods: 


Ike A quasi-steady procedure developed 
by McCarthy at David W. Taylor 
Naval Ship Research and Development 
Center (DTNSRDC), 

Be An unsteady lifting surface theory 
developed by Tsakonas and his 
colleagues at Davidson Laboratory, 

336 An unsteady lifting surface theory 
developed by Kerwin and Lee at MIT, 

4. A refinement by Kerwin to the 
method of Kerwin and Lee to con- 
sider the inclination of the pro- 
peller slipstream. 


In inclined flow, all four of the 
calculation methods evaluated consis- 
tently underpredicted the experimental 
values of the periodic single-blade 
loads. The method of Kerwin, which 
considers the inclination of the slip- 
stream relative to the propeller axis, 
produced the best correlation with ex- 
perimental values. These correlations 
show that the inclination of the pro- 
peller slipstream relative to the pro- 
peller axis significantly influence the 
periodic single-blade loads. The impor- 
tance of this inclination increases with 
increasing timeaverage loading. 


In longitudinal flow, all of the 
calculation procedures predicted period- 
ic singleblade axial force at design 
advance coefficient to within 20 percent 
of the experimental values, but agree- 
ment waS not as good at off-design ad- 
vance coefficients. The method of 
Kerwin and Lee produced the best overall 
correlation with experimental results 
considering both amplitude and phase of 
the periodic single-blade axial force 
over a range of advance coefficients. 


NOMENCLATURE 
R 

Ar Expanded area, Z| cdr 
Th 

Ao Propeller disk area, 


mD2/4 


c Blade section chord length 


Fx,y,2 


V 
VA 


Vr (©, 8y) 


Ve (Ll, Ow) 


Vx (rl, Ow) 


VvM 


Propeller diameter 


Force components on blade 
in x,y,z directions 


Advance coefficient, 
J=Va/nD 


Torque coefficient, 
Q/ (on2p9) 


Thrust coefficient, 
T/ (pn2p4) 


Moment components about 
x,y,z axes loading on one 
blade 


Propeller revolutions per 
unit time 


Propeller blade section 
pitch 


Time-average propeller 
torque arising from 
loading on all blades, 
-2My 


Radius of propeller 


Radial coordinate from 
propeller axis 


Time-average thrust of 
propeller, positive 
forward, ZF x 


Model speed 


Propeller speed of advance 


Radial component of wake 
velocity, positive towards 
hub 


Tangential component of 
wake velocity, positive 
clockwise looking 
upstream for left-hand 
propeller, positive 
counter-clockwise looking 
upstream for right-hand 
propeller 


Longitudinal component of 
wake velocity, positive 
forward 


Volume mean longitudinal 
velocity through propeller 
disk determined from wake 
survey 


(Op mM)n 


Subscripts: 
d 


EXP 


Wake fraction determined 
from volume mean 
longitudinal velocity 
through propeller disk 
determined from a wake 
survey, (V-Vym)/V 


Coordinate axes 
Number of blades 


Angular coordinate used 
to define location of 
blade and variation of 
loads, from vertical 
upward, positive 
counterclockwise looking 
upstream for left-hand 
propeller, positive 
clockwise looking 
upstream for right-hand 
propeller, 0= -6y 


Angular coordinate of 
wake velocity, from 
upward vertical, positive 
clockwise looking up- 
stream for left-hand 
propeller, positive 
counterclockwise looking 
upstream for right-hand 
propeller, ®w =-6 


nth harmonic force 
coefficient, 


(F) n/ (onvaD3) 


nth harmonic moment 
coefficient, 


(M) n/ ( onVaD4) 
Mass density of water 


nth harmonic phase angles 
of F,M based on a cosine 
series, 
WAG A N 
(F,M) = (F,M)+ © (F,M)n. 
n=1 
cos {né-(¢F,m) n} 


Inclination of propeller 
shaft 


Value at the near design 
condition 


Value of experimental 
measurements 


Value of hub radius 


Value of nth harmonic 


REF Value at specified 
reference advance 
coefficient 

X,VY,Z Component in x,y,z 
direction 


Superscripts: 


- Time-average value per 
revolution 


INTRODUCTION 


The once-per-revolution (or shaft 
frequency) propeller blade loads are 
important for consideration of propeller 
blade fatigue strength (1,2,3,4). 
Therefore, an investigation was 
undertaken at DTNSRDC to evaluate the 
periodic blade loads on propellers in 
inclined flow and in circumferentially 
non-uniform longitudinal velocity fields. 


Systematic experiments were 
conducted on three model propellers (see 
Table I) operating in inclined flows 
generated by shaft inclinations of 10, 
20, and 30 degrees, and in a once- 
per-revolution circumferentially 
non-uniform longitudinal velocity field 
produced by an upstream wire-grid 
screen. The inclined flow produces a 
once-per-revolution variation of the 
tangential component of velocity in the 
propeller plane; therefore, the inclined 
flows are sometimes referred to as 
tangential wake patterns. 


TABLE I — CHARACTERISTICS OF PROPELLERS 
4661, 4710, AND 4402 


PROPELLER NUMBER | 4661 | 4710 | 4402 | 


Diameter, D (feet) 


Rotation 


Number of Blades, Z 


Hub-Diameter Ratio, D,,/D 


Expanded Area Ratio 
Design Advance Coefficient 


Pitch-Diameter Ratio (P/D) 
atr/R =0.7 


The experimental results were 
correlated with predictions by the 
following theoretical methods: 


Ib A simple quasi-steady procedures 
developed by McCarthy (5) at 
DTNSRDC which utilizes the open 


BLADE ON WHICH 
FORCES ARE MEASURED 


AHEAD (+v) 


PROPELLER 
CENTERLINE 


COORDINATE SYSTEM FOR RIGHT HAND PROPELLER 
(PROPELLERS 4710 and 4402) 


THE COORDINATE SYSTEM 
ROTATES WITH THE PROPELLER 


F,,2 BLADE ON WHICH 
FORCES ARE MEASURED 


PROPELLER 
CENTERLINE 


COORDINATE SYSTEM FOR LEFT HAND PROPELLER 
(PROPELLER 4661) 


Fig. 1 — Components of Blade Loading 


water characteristics of the pro- 
peller, 

Ao A linearized lightly-loaded 
unsteady lifting surface theory 
developed by Tsakonas, et al (6,7) 
at Davidson Laboratory, 

3. A moderately-loaded unsteady 
lifting surface theory developed by 
Kerwin and Lee (8) at MIT, 

ay A refinement by Kerwin (9) to the 
method of Kerwin and Lee which 
considers the inclination of the 
propeller slipstream for opera- 
tion in inclined flow. 


BACKGROUND 


In previous work unsteady loads 
(see Figure 1) were measured on a single 
blade of model propellers operating 
behind model hulls with open-shaft 
transom sterns (1,2,3,4). It was found 
that for these vessels the unsteady 
blade loads were produced primarily by 
the tangential component of the wake 
velocity; see Figure 2. For the steady 
ahead conditions, the experimentally 
obtained unsteady loads were correlated 
with unsteady loads deduced from strains 
Measured on the respective full scale 
propeller blades and with predictions 
based on unsteady lifting surface theory 
developed by Tsakonas, et al (6,7), and 
the quasi-steady method of McCarthy (5). 


For each of these propeller-hull 
combinations, the circumferential 


variation of loading determined from the 
model experiments agreed fairly well 
with full-scale data, but was substan- 
tially larger than the predictions by 
the two theoretical methods. 


In an attempt to isolate the 
reasons for the large discrepancy 
between these experimental and 
theoretical results, the following 
experimental work on model propellers 
was undertaken: 


V, Positive Downstream 


(Looking Upstream) 


Fig. 2 — Sign Convention for Ship Wake 
Velocity Components 


Ike Measurement of the six components 
of blade loading under various 
idealized conditions including: 


a. In inclined flow without an 
upstream hull. 
b. In a once-per-revolution 


variation of longitudinal flow 
produced by an upstream 
wire-grid screen. 

Cre Over a range of propeller 
advance coefficients. 

d. Over a range of clearances 
between the propeller tip and 
a flat plate above the 
propeller. 

2. Measurement of the field point 
velocities induced by a propeller 
operating in inclined flow. 

Bie Measurement of the pressure 
distribution on the surfaces of the 
blades of two propellers operating 
in inclined flow. 


This paper presents the results of 
the blade loading experiments under 
idealized conditions, and correlation 
with available theories including one 
which considers the influence of 
slipstream for operation in inclined 
flow. These results supercede the 
preliminary results presented in 
Reference 10. More details of these 
results will be presented in future 
DTNSRDC Reports. The field point 
velocity measurements were reported in 
Reference 11. The results of the blade 
surface pressure measurements will be 
presented in a future DTNSRDC Report. 


EXPERIMENTAL TECHNIQUES 


The experiments were conducted on 
DTNSRDC Carriage 2. The propellers were 
mounted at the front of the downstream 
drive system containing the dynamometer 
and the drive motor. This system was 
the same as that used in the experi- 
ments described in detail in References 
2 and 3 except that the downstream body 
was modified so that it could be 
operated fully submerged. The 
modifications included a waterproof 
housing for the drive motor, waterproof 
electrical cables and connectors, 
removal of the upper apron which had 
extended the sides of the boat, and the 
addition of a non-waterproof top to the 
boat. 


The sensing elements were flexures 
to which bonded semi-conductor 
strain-gage bridges were attached. 

Three flexures were necessary tO measure 
all six components of forces and 
Moments. Flexure 1 measured Fy and 

My, Flexure 2 measured Fy and My, 

and Flexure 3 measured Fz and Moi 


see Figure 1. The flexures were mounted 
\ 


inside a propeller hub specifically 
designed for these experiments. Only 
one flexure could be mounted at a time 
because of space limitations, and this 
necessitated three duplicate runs for 
each condition. The flexure calibration 
procedure was identical to that 
described in References 2 and 3. 


Wake Screen and Wake Survey 


A 16-inch (0.405 m) diameter wake 
screen designed to produce a dominant 
first harmonic of the circumferential 
variation of the axial velocity, was 
used in the experiments. This screen 
consisted of simply a base screen with 
one overlay screen and no support 
members in the region through which flow 
is drawn into the propeller disk. 


Wakes were measured for three 
different flow regimes: 
Ls Shaft inclination = 10 deg with 
no wake screen, 
De y 20 deg with no wake screen, 
3a w 0 deg with the wake screen. 


Figure 3 shows the velocity 
component ratios at the radial position 
near r=0.7R at which the various wake 
surveys were conducted. The 
circumferential variations of the 
velocity component ratios in the axial, 
tangential, and radial directions 
contained primarily first harmonic 
components. 
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Fig. 3 — Distribution of Wakes in 
Propeller Disks 
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Fig. 3 — (Continued) 


Figure 3 shows that the variation 
of the longitudinal velocity component 
ratio Vx/V behind the wake screen is 
approximately a step function, as was 
expected. However, the high velocity 
region covers a smaller portion of the 
propeller disk than the low velocity 
region, whereas it was anticipated that 
the two velocity regions would each 
cover essentially half of the disk. 
This suggests that there is some cross 
flow between the propeller and the 
screen. This non-symmetry of the wake 
profile causes no problem regarding the 
objectives of the present investigation 
since it remains the same for all 
experimental conditions. 


Experimental Conditions and Procedures 


The propeller was located 
approximately 16 in (0.41 m) downstream 
of the screen for all conditions for 
which the screen was used. A 4.0 ft by 
2.0 ft by 0.375 in (1.22 m by 0.61 m by 
9.5 mm) plate was positioned parallel to 
the flow above the propeller for all 
experimental conditions. For Propeller 
4402 at 10 degrees inclination, the 
plate clearance from the tip of the 
propeller was varied from 1.0 in (0.025 
m) to 20 in (0.51 m). For all other 
conditions the plate clearance was 
approximately 20 in (0.51 m). With this 
clearance, the plate was 2 in (0.051 m) 
below the water surface. The plate at 
this location was used to minimize any 
possible effects of the free surface on 
the propeller performance. 


The experimental conditions are as 
follows: 
ave No Wake Screen, w = 10 deg, 

Propellers 4661, 4710, and 4402. 
Qi No Wake Screen, = 20 deg, 

Propellers 4661, and 4710. 

3% Wake Screen, » = 0 deg, 

Propellers 4661 and 4710. 

4. No Wake Screen, wW = 30 deg, 

Propeller 4661. 

For each condition, the carriage 
speed was held constant and propeller 
rotation speed, n, was varied to obtain 
the desired range of advance 
coefficient, J. For each condition, the 
carriage speed for the blade loading 
Measurements was the same as the 
carriage speed for the corresponding 
wake survey. 


The first three conditions 
represent the more important types of 
experiments performed. The wake surveys 
were conducted at these conditions. The 
fourth condition was run to obtain 
limited data at a very large value 
of ~ so that the variation of loads 
with ~ could be better defined. Due to 
limited time, however, only the Fy and 
My components were measured on Pro- 
peller 4661, and no measurements were 
made on Propeller 4710. In addition, no 
wake survey was conducted at this 
condition due to limitations of the wake 
Survey apparatus. 


For each of these conditions, 
experiments were conducted over a range 
of advance coefficient J. For each 


condition, between 20 and 40 values of 
advance coefficient were evaluated in 
order to well define the manner in which 
the loading varies with advance 
coefficient. In general, each 
experimental condition was run three 
times, once for each of the three 
flexures. This was necessary in crder 
to obtain all six components of force 
and moment. Data were also collected 
with the propellers rotating in air at 
0, 10, 20, and 30 degrees of shaft 


inclination over the range of rotational 
speeds evaluated in water. This 
provided the gravitational and centri- 
fugal loads that were later subtracted 
from the total propeller blade loads 
measured in water. 


Data Acquisition and Analysis 


Data were collected and analyzed in 
the same manner as in the steady ahead 
runs described in References 1 through 
4. For each flexure, an Interdata 70 
Computer was used to collect and average 
the force and moment data for each 
4-degree increment of blade angular 
position over 200 to 300 propeller 
revolutions. Propeller rotation speed, 
n, and model (or carriage) velocity, V, 
were recorded for every set of two 
propeller revolutions, and averaged over 
the period during which data were 
collected. 


At the end of each run, the 
computer analyzed and printed the data. 
The average force and moment signals for 
each 4-degree angular position were 
printed along with the average model 
velocity and propeller rotation speed 
for the run. The standard deviation of 
the accumulated data for the run was 
also calculated for V, n, and the force 
and moment signals at each position. A 
harmonic analysis was performed on the 
force and moment data providing the mean 
signal and amplitude and phase of the 
first 16 harmonics. 


First, the experimental agenda was 
completed for all three flexures at each 
Major condition (each propeller, wake 
screen, and value of ~) except wp = 30 
deg. Then the corrections for 
interactions were calculated at five 
values of J at which measurements with 
all three flexures were made using the 
same procedures outlined in Reference 1 
through 4. After the correction for 
interactions, the centrifugal and 
gravitational loads were subtracted from 
the total experimental loads to obtain 
the hydrodynamic loads. 


Accurac 


The accuracy of the experiment was 
generally similar to that of the steady 
ahead runs in References 1 through 4. 
During a single run, where signals were 
averaged over many revolutions, the 
standard deviation of the measured 
forces and moments maintained a similar 
error (95 percent confidence band) of +5 
to +10 percent of the averaged signal at 
each angular position. For a given run 
the average error in model speed V was 
approximately +1.0 percent, and the 
error in rotational speed n was less 
than +0.01 percent 


Except for the fluctuations in 
signals occurring in a given run, the 
accuracy of the data is indicated by the 
repeatability of an identical condition. 
An effort was made to set experimental 
conditions identically on repeat runs; 
however, some variation was unavoidable 
due to manual setting of the model speed 
and propeller rotation speed, random 
zero shift in data channels, and minor 
variations in basin water conditions. 
This resulting error is approximately +2 
percent of the measured first harmonic 
loading component. A more detailed 
error analysis based on the repeata- 
bility of the hydrodynamic, centrifugal 
and gravitational loads will be presented 
in a future DTNSRDC Report. 


EXPERIMENTAL RESULTS 


Loading Components 


The basic loading components are 
shown in Figure 1. For a right-hand 
propeller the sign convention follows 
the conventional right-hand rule with 
right-hand Cartesian coordinate system. 
For a left-hand propeller all the loads 
are the same, but for this case the sign 
convention follows a left-hand rule with 
a left-hand Cartesian coordinate system. 


Each component of loading is 
generally represented as a variation of 
the instantaneous values with blade 
angular position, 6, and as a Fourier 
series in blade angular position in the 
following form: 


epee N 
F,M(0) = (F,M)+ > (F,M), cos{nd-(d¢ m),} (2) 


n=1 


In general, the loads consist of 
hydrodynamic, centrifugal and gravita- 
tional components. However, in this 
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Fig. 4 — Experimental Variation of Loading Components with Blade Angular Position 


paper only the hydrodynamic component of 
blade loading is presented. Centrifugal 
and gravitational loads were measured 
soley to permit the hydrodynamic loads 
to be determined by subtracting the 
centrifugal and gravitational loads from 
the total experimental loads. 


Centrifugal and Gravitational Loads 


Centrifugal and gravitational loads 
were determined from air-spin experi- 
ments with each flexure over a range of 
rotational speed n for ¥ = 0, 10, 20, 
and 30 deg. The centrifugal load, which 
is a time-average load in a coordinate 
system rotating with the propeller, 
should vary as n2 and be independent 
of ~. The time-average experimental 
data followed these trends closely. The 
gravitational load, which is a first 
harmonic load in a coordinate system 
rotating with the propeller, should vary 
with and be independent of n. The 
first harmonic experimental data 
followed these trends closely. 


The centrifugal and gravitational 
loads measured during these experiments 
agreed with the values previously 
reported in References 3 and 4 and the 
gravitational loads agreed with values 
deduced from the weights of the blades 
and associated flexures. Therefore, 
these results will not be repeated here. 


Variation of Loads with Angular Position 


Figure 4 shows the variations of 
the loading components with blade 
angular position for Propeller 4661 for 
operation with 10 degrees shaft 
inclination with no screen and for 
operation behind the wake screen with no 
shaft inclination. These plots present 
typical results obtained in inclined 
flow and behind the wake screen. For 
other conditions evaluated in inclined 
flow and behind the wake screen the 
trends are basically the same as shown 
in Figure 4, but the magnitude of the 


unsteady loading varies with 
experimental conditions. 


The variations of the primary load 
components, Fy and M, with blade 
angular position in inclined flow and 
behind the wake screen follow patterns 
approximately similar to the tangential 
and the longitudinal velocity profiles 
of the respective wakes shown in Figure 
3. The variation of the loads behind 
the wake screen implies that the high 
velocity region covers a smaller portion 
of the propeller disk than the low 
velocity region, which agrees with the 
measured wake data. 


Effect of Plate Clearance 


Figure 5 shows that the periodic 
blade loads in inclined flow are fairly 
insensitive to the presence of a nearby 
flat boundary parallel to the flow for 
tip clearance to propeller diameter (tip 
clearance ratio) as small as 0.1. The 
ratios of the periodic to the time- 
average loading components at tip 
clearance ratio of 0.1 are no greater 
than 10 percent larger than the 
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Fig. 6 — Experimental Variation of First Harmonic Blade Loads 
On Propeller 4661 with Advance Coefficient 


corresponding values at very large tip 
clearance ratios for Propeller 4402 with 
an inclination angle ~ of 10 degrees. 


Effect of Advance Coefficient 


Figure 6 presents the trends of the 
variations of the first harmonics of the 
Bree load components Fy, M 

and My with J for Propeller 4661 
aye the three primary experimental 
conditions described in the section on 
Experimental Conditions and Procedures. 
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In Figure 6 the force and moment 
components are nondimensionalized as 
follows: 


(Eee) 


ag x,y'1 
Meaehn @ = (2) 
pnV,D 
a af (My ys (3) 
ESD Sar 
pnV,D 
where the subscript 1 represents the 
first harmonic component, and Va = 
V(l-wym) is the speed of advance based 
on the volume mean wake. This form of 
nondimensionalization was used in part 
to verify Wereldsma's argument (12) that 


for a given propeller in a given wake 
pattern the circumferential variation of 
the hydrodynamic loading varies 
approximately as nVa; i.e., for a 

given value of nVaq the circum- 
ferential variation of hydrodynamic 
loading ts insensitive to J. This form 
of nondimensionalization is different 
only by a factor of the advance coeffi- 
cient, Ja, from the conventional form 
of nondimensionalization; i.e., 
Fx,y/en2p4 and My, y/en? 2p5. 

The conericrent: (EX, y)l 

and (mx,y)1 in Figure 4s are all 
normalize by the respective 
coefficients at design J to illustrate 
the relative sensitivity of the loading 
coefficients to J, for different wakes. 


The results presented in Figure 6 
show that in tangential flow ‘(px,y)1 
and (mx,y)1 generally decreases with 


increasing J. K(My)1 is the most 
sensitive to J, and (py), is the 

least sensitive to J. Bxcept for the 

Fy component, these data do not 

closely follow the trends indicated by 
Wereldsma (12). For the longitudinal 
wakes behind the screen, the slopes 

of «(pF m) with increasing J are less 
negative than in inclined flow. The 
relative slopes among the four 
components are generally the same for 
each propeller-wake combination, except 
for some cases at greater than design J. 


Figure 7 compares the variations 
of (px)1 with J obtained in the 
present experiment in inclined flow with 
variations obtained in previous exper- 
iments in which the unsteady blade loads 
were determined over a range of J in 
inclined flow (13,14,15,16). Each of 
these sets of experimental results show 
that (py), decreases substantially 
with increasing J. 


Figure 8 compares the variation 
in k(Fx)1 with J obtained in the 
present experiment behind the wake 
screen with the blade frequency K(Fx)n 
obtained in previous experiments in 
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Fig. 7 — Comparison of First Harmonic Blade 
Loading Coefficients in Tangential Wakes 


which the unsteady bearing forces were 
measured in longitudinal wakes produced 
by screens in a closed-jet water tunnel 
(17,18,19). The results in References 
17, 18, and 19 are for the blade fre- 
quency component of ‘py, whereas in 

the present investigation the shaft 
frequency component of Kkpy was 

measured. However, the mechanism for 
generating the unsteady blade loading in 
axial flow appears to be independent of 
the harmonic of shaft rotation. Each of 
the three previous sets of experiments 
shows that k(px)n decreases with 
increasing J to some minimum value near 
design J, and then increases with 
further increase in J. Figure 8 
illustrates that variations of k(py)1 
with J obtained in the present 
experiment generally follow the trend of 
the data in References 17, 18 and 19. 

In general, “(Px)n is somewhat less 
sensitive to in longitudinal wake 
patterns than it is in tangential wake 
patterns. 
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CORRELATION BETWEEN EXPERIMENTAL RESULTS 
AND THEORETICAL PREDICTIONS 


Theoretical Methods 


The experimental results were 
correlated with predictions based on the 
following methods: 


IG The quasi-steady method developed 
at DTNSRDC by McCarthy (5); 
Computer Program QUASI. 

aie The procedure developed 
Laboratory by Tsakonas, 
based on lightly-loaded 
lifting surface theory; 
Program PPEXACT. 

3}6 The procedure developed at MIT by 
Kerwin and Lee (8) based on 
moderately-loaded unsteady lifting 
surface theory; Computer Program 
PUF2. 

4. A refinement of the method of 
Kerwin and Lee (8) developed at MIT 
by Kerwin (9), to allow the axis of 
the propeller slipstream to depart 
from the propeller axis for opera- 
tion in inclined flow; Computer 
Program PUF2IS. 


at Davidson 
et al (6,7) 
unsteady 
Computer 


The procedure developed by McCarthy 
(5) is a simple quasi-steady procedure 
utilizing the open water characteristics 
of the propeller. It is assumed that 
the thrust and torque developed by the 
propeller blade at any angular position 
in a circumferentially nonuniform wake 
is the same as would be produced by the 
propeller blade if it were operating 
continuously at the advance coefficient 
J and rotational speed n based on the 
local wake at the angular position 
corresponding to the mid-chord of the 70 


percent radius. It is further assumed 
chat the instantaneous thrust and torque 
can be adequately estimated by entering 
the propeller open water characteristics 
at the values of J and n based ona 
weighted average over the propeller 
radius of the wake at the local blade 
angular position. This simple method 
can be expected to yield reasonable 
results only if the reduced frequency of 
interest is low, the propeller projected 
skew is small relative to the wave 
length of the pertinent wake harmonic, 
and the wake harmonic corresponding to 
the force harmonics of interest does not 
vary substantially in amplitude or phase 
with radius. These conditions are met 
for propellers and wakes being evaluated 
in the present paper. 


The procedure developed by 
Tsakonas, et al (6,7), is based on the 
linearized unsteady lifting surface 
theory for a lightly-loaded propeller 
using an acceleration potential. The 


10 


numerical procedure applies the mode 
approach and collocation method in 
conjunction with the "generalized lift 
operator" technique. This procedure 
based on the frequency domain assumes 
that a given harmonic of blade loading 
depends upon the corresponding harmonic 
of the wake velocity normal to the blade 
chord line, independent of whether the 
normal velocity results from axial or 
tangential components of the wake. 

The blade load does not depend upon the 
radial variation of the circumferential 
mean wake. The shed and trailing 
vortices are assumed to lie on an 
"exact" helicoidal surface of constant 
pitch extending to downstream infinity 
determined by the propeller rotational 
speed and a single axial inflow 
velocity. The axis of this helicoidal 
surface coincides with the propeller 
axis, regardless of the inclination of 
the propeller shaft to the incoming 
flow. This method does not consider the 
contraction or roll-up of the propeller 
slipstream. All geometric characteris-— 
tics of the propeller are considered, 
except rake, camber, and thickness which 
are assumed to be zero. 


Valentine (20) developed a 


refinement to the method of Tsakonas, et 


al (6,7) for operation in slightly 
inclined flow. This refinement, in 
effect, replaces the "exact" helicoidal 


wake whose axis coincides with the 
propeller axis with a slightly distorted 
one in the direction of the inflow 
velocity. The refinement, which is 
incorporated as a perturbation for small 
inclination angles, relates the unknown 
loading at the first harmonic of shaft 
frequency with the loadings at the mean 
and second harmonics of shaft frequency 
evaluated without the distorted helicoi- 
dal wake. All other assumptions of the 
Method of Tsakonas et al, are retained. 


Calculations made by Valentine (20) 
showed that his modification to the 
method of Tsakonas et al, did not 
significantly improve the poor 
correlation between this method and 
experimental blade loads in inclined 
flow. Further, with the modifications 
by Valentine the disagreement in phase 
between predicted and experimental blade 
loads in inclined flow were even worse 
than the agreement obtained with the 
method of Tsakonas et al. Valentine 
concluded that the effects of shaft 
inclination required a moderately-loaded 
propeller theory, and could not be 
adequately calculated based on the 
lightly-loaded formulation of Tsakonas 
et al. Therefore, the theory of 
Tsakonas et al, as modified by Valentine 
was not correlated with the experimental 
results in the present paper. 


The numerical method developed by 
Kerwin and Lee (8) is based on a 
linearized lifting-surface theory in the 
time domain. The propeller blades are 
represented by a spanwise and chord- 
wise distribution of discrete line 
vortex and source elements located on 
the exact camber surface of the blade. 
Thus the geometric complications of 
skew, rake and radial variation in pitch 
are readily accommodated. The trailing 
vortex wake is permitted to contract and 
roll up, and the effect of vortex sheet 
separtation from the blade tip is taken 
into consideration. The inflow velocity 
to the propeller may have radially and 
circumferentially varying axial, 
tangential, and radial components, and 
may therefore give rise to both steady 
and unsteady blade loading. This 
method, like the method of Tsakonas et 
al, assumes that the axis of the 
propeller slipstream coincides with the 
propeller axis. 


Kerwin (9) developed a refinement 
to the method of Kerwin and Lee (8) for 
operation in inclined flow where the 
slipstream is not axisymmetric about the 
propeller shaft. This refinement 
entails a more realistic representation 
of the path of the propeller 
slipstream. In this method the axis of 
the slipstream coincides with the 
propeller axis immediately behind the 
propeller, and coincides with the 
direction of the mean inflow far 
downstream in the ultimate wake. A 
simple function is assumed for the 
slipstream axis in the wake region 
between the propeller and the ultimate 
wake. Due to the asymmetry, the 
position of the trailing vortex wake 
relative to a blade oscillates with a 
once-per-revolution fundamental fre- 
quency, thus giving rise to unsteady 
induced velocities normal to the blade 
surface and thereby unsteady blade 
loadings of the same frequency. The 
strength of the vorticity in the wake, 


and thus the induced once-per-revolution 
variation of loading on the blades, is 
dependent upon the time-average load- 
ing of the propeller. All other 
characteristics of the method of Kerwin 
and Lee except for roll-up are retained, 
including the flexibility for the 
trailing vortex wake to contract, and 
allowance for the effect of vortex sheet 
separation from the blade tip. 


Correlations in Tangential Wakes 


Figures 9 and 10 present the 
amplitudes and phases, respectively, of 
the first harmonic loads on Propellers 
4661, 4710, and 4402 operating in the 
various tangential wakes over a range of 


JEL 


advance coefficient J. These figures 
present both experimental results and 
predictions based on the various 
theoretical methods described in the 
preceding section. For Propeller 4661 
with 10 degrees shaft inclination, the 
first harmonic components of Fy, My, 
Fy, and My are presented. For other 
conditions only the Fy component is 
shown. The results in Figure 9 and 10 
are summarized in Table II for design J. 


In the tangential wakes, a 
consistent variation from the 
experimental data occurred in the 
theoretical predictions for the three 
propellers evaluated. With few isolated 
exceptions, all of the theoretical 
methods underpredicted all of the 
loading components throughout the range 
of conditions evaluated. In general, 
the correlations are better for the Fy 
and My components than for the Fy 
and My components. The Fy and M 
components are the more important 
components since, in general, these are 
larger than the Fy and My components. 


The quasi-steady method of McCarthy 
(5) underpredicted the amplitude of 
(Fy)1 by approximately 10 to 30 
percent of the experimen tal values with 
closer agreement at higher values of J. 
The quasi-steady predictions main- 


tained a similar trend to the 
experimental data for Fy and My 
components for all conditions in 
tangential wakes; i.e., both the 
experimental and predicted values of 
K(Px) 1 and (Mx) decrease with 
increasing J. or a given propeller the 
predicted slope of the K(Px)1 and 
K(Mx) 1 Curves (not shown) with J 
increases with increasing amplitude of 
the tangential wake. No predictions of 
(Fy)1 and (My)1] were made with 

the quasi-steady procedure since this 
method does not predict the radial 
center of the load. These components 
could be predicted by the quasi-steady 


method by assuming a radial position of 
application of the load. 


The phases presented in Figure 10 
show that the quasi-steady method of 
McCarthy predicts that the maximum 
values of the loading components will 
occur at approximately 10 to 20 degrees 
of blade angular position before the 
experimentally determined angles. 


Predictions by the unsteady theory 
of Tsakonas et al (6,7) did not agree 
well with experimental results in 
tangential flow. This theory predicts 
that the amplitudes of all four loading 
coefficients increase with increasing J, 
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which is opposite to the trend of the 
experimental data. Since this is a 
lightly-loaded linear theory, it was 
anticipated that the predicted unsteady 
loading coefficients «(p my1, would be 
essentially independent Off JAN 
earlier version of the unsteady lifting 
surface method developed by Tsakonas and 
his associates (21) at Davidson 
Laboratory, in which the propeller 
helical wake was approximated in a 
staircase manner did predict that the 
periodic loading coefficients ‘(F,m) 

in axial wakes were essentially 
independent of J (17). It is not clear 
why the method of Tsakonas et al (6,7) 
predicts that kp m)1 increases 
substantially with increasing J. 


In general, the predictions by the 
method of Tsakonas et al are closer to 
the experimental results at high advance 
coefficients than they are at low 
advance coefficients. The predicted 


amplitudes are approximately 65 percent 
of the experimental value at design J 
and 30 to 35 percent of the experimental 
values at the lowest values of J 
evaluated. The finding that this method 
predicts periodic loads in inclined flow 
that are significantly less than the 
experimental values is consistent with 
the results presented in References 1 
through 4. 
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The unsteady theory of Tsakonas et 
al predicts that the maximum values of 
the loading components will occur at 


approximately 15 to 20 degrees of blade 
angular position before the 
experimentally determined angles; see 


Figure 10. 


The unsteady theory of Kerwin and 
Lee (8) showed, in general, somewhat 
better agreement with experimental data 
over a range of advance coefficient than 
either the quasi-steady method of 
McCarthy or the unsteady method of 
Tsakonas et al. This applies to both 
amplitudes and phases. The trends of 
the predictions of the amplitudes over a 
range of advance coefficients were very 
similiar to the trends of the predic- 
tions by the method of McCarthy. The 
method of Kerwin and Lee underpredicted 
the periodic loads, by approximately 5 
to 20 percent of the experimental 
values, with closer agreement being 
obtained for the Fx and My 
components at the higher values of J. 


The unsteady theory of Kerwin and 
Lee predicts the phases to within 5 to 
15 degrees of the experimental values 
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TABLE II — COMPARISON OF EXPERIMENT AND THEORY IN INCLINED FLOW — AMPLITUDES AND PHASES 
OF FIRST HARMONIC LOADING COEFFICIENTS AT DESIGN J 


PROPELLER 4661, Y = 10 DEG, J= 1.14 
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for most conditions evaluated. In 
general, the agreement is better for the 
higher values of advance coefficient J; 
i.e., J near or higher than the design 
value. 


The method of Kerwin (9), which is 
a refinement to the method of Kerwin and 
Lee to account for the effect of the 
inclined slipstream, showed substantial 
improvement over the predictions of 
amplitudes and phases by the method of 
Kerwin and Lee at values of J which are 
substantially less than the design J. 
For higher J values, i.e., J near or 
higher than the design value, where the 
method of Kerwin and Lee was in closer 
agreement with experimental results, the 
improvement by the method of Kerwin was 
smaller. This implies that the 
influence of the non-axi-symmetric wake 
due to shaft inclination on the 
fluctuating forces becomes larger at 
lower values of advance coefficient. 
This trend is expected, as pointed out 
by Kerwin (9), since the effect of wake 
asymmetry is the combined result of the 
time-average loading and the oscillating 
position of the trailing vortex wake 
relative to a point fixed on the blade. 
At low advance coefficients, the mean 
loading is higher, and the pitch of the 
trailing vortex wake is smaller. This 
smaller pitch of the propeller wake 
tends to increase the influence of the 
fluctuating wake position relative to 
the blade on induced velocities since 
the downstream wake is closer to the 
propeller blades. 


In summary, these correlations show 
that the inclination of the propeller 
slipstream relative to the propeller 
axis can significantly influence the 
periodic loads on the propeller blades. 
The importance of this inclination 
increases with increasing time-average 
loading. Of the four methods evaluated, 
the method of Kerwin, which accounts for 
the inclination of the slipstream, gives 
the best predictions of the amplitudes 
and phases of the periodic blade loads. 
The method of Tsakonas et al gives the 
worst prediction of the amplitudes of 
the periodic loads. Further, this 
method failed to predict the trend of 
the variation of the loading 
coefficients with advance coefficient. 


Figure 1]1 presents the variation of 
the periodic loads with inclination 
angle for Propeller 4661 near design J. 
These results show that the rate of 
increase of the periodic loads with 
inclination angle increases with 
increasing inclination angle. Figure 11 
also shows that the correlation between 
experimental resulsts and analytical 
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method of Kerwin, which considers the 
inclination of the propeller slipstream, 
is essentially independent of the angle 
of inclination up to 30 degrees. 


PROPELLER 4661 
Oo——O © Experiment 


&r-—-- © Theory (Kerwin-PUF2IS) 


J = 1.062 


INCLINATION ANGLE, y (DEG) 


Fig. 11 — Variation of Periodic Blade Loads with Shaft Inclination 


Correlations in Longitudinal Wakes 


Figures 12 and 13 present the 
amplitudes and phases, respectively, of 
the first harmonic of the Fy component 
on Propellers 4661 and 4710 in the 
longitudinal wake patterns generated by 
the upstream wake screen over a range of 
advance coefficient J. The 
circumferential variation of the 
longitudinal component of the incoming 
velocity is essentially a once-per- 
revolution variation with low velocity 
on the left half of the propeller disk 
and high velocity on the right side of 
the propeller disk looking upstream; see 
Figure 3. Figures 12 and 13 also 
present the correlation between experi- 
mental results and predictions based on 
the three theoretical methods described 
previously.* The results in Figures 12 
and 13 are summarized in Table III for 
design J. 


*For 
the 
the 
Lee 


longitudinal inflow where yp = 0, 
method of Kerwin (9) is identically 
same as the method of Kerwin and 
(8) . 
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Fig. 12 — Amplitudes of First Harmonic Blade Loads in Longitudinal Wakes; Correlation 
Between Experiment and Theory 
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Fig. 13 — Phases of First Harmonic Blade Loads in Longitudinal Wakes; Correlation 
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TABLE III - COMPARISON OF EXPERIMENT AND THEORY IN 
LONGITUDINAL FLOW — AMPLITUDES AND PHASES 
OF FIRST HARMONIC LOADING COEFFICIENTS 
AT DESIGN J 


Kies 


Kex)4 


K (pdt, EXP K(Ex)1, EXP 


amPLITUDE | 1000- 
(Fx)1 


Experiment 


McCarthy 
(QUASI!) 


Tsakonas et al 
(PPEXACT) 


Kerwin and Lee 
(PUF2) 


Experiment 


McCarthy 
(QUASI) 


Tsakonas et al 
(PPEXACT) 


Kerwin and Lee 
(PUF2) 


The experimental curve presenting 
“(Px)1 versus advance coefficient J 
shows a concave shape with a minimum 


value near design J; see Figure 12. The 
same trend was observed for other loading 
components (not shown). All of the theo- 
retical methods are in good agreement 
with experimental (py) ] near design 

J; however, they did not, in general, 
agree well with experimental (py) 1 


over a wide range of J. The correla- 
tions between experiments and theore- 
tical predictions are somewhat different 
for the two propellers. 


The predictions by the quasi-steady 
method of McCarthy showed different 
trends for different propellers. For 
Propeller 4661, the predicted values of 
and *(px)1 and 90k/dJ increase with 
increasing J. For Propeller 4710, how- 
ever, the predicted values of X(py)1 
and 3k/dJ decrease with increasing J. 
The trends of the predictions by the 
quasi-steady method of McCarthy are quite 
sensitive to small changes in the slopes 
of the open water curves; i.e. Km 
versus J and Kg versus J. 


The phases shown in Figure 13 show 
that the quasi-steady method of McCarthy 
predicts that the maximum value of 
(Fy)1, will occur at approximately 


PROPELLER 4661 | PROPELLER 4710 
CONDITION vy =0 DEG ) =0 DEG 
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25 and 10 degrees of blade angular 
position after the experimentally 
determined angles for Propellers 4661 
and 4710, respectively*. 


The unsteady theory of Tsakonas et 
al (6,7) predicts a similar trend of 
K(Fx)] with J in axial and tangential 
wakes; i.e., the predicted (py) 1 
tends to increase with increasing J 
and 0k/dJ to decrease with increasing 
J. These predicted trends are similar 
because this method considers only the 
component of wake resolved normal to the 
blade pitch line, for either longi- 
tudinal or tangential wakes. Thus, as 
discussed previously, this method does 
not distinguish between longitudinal and 
tangential wakes. 


The trend of K(py)] with J 
predicted by the method of Tsakonas et 
al does not agree with the experimental 
results. At low values of J; i.e., high 
time-average loading coefficient, this 
method gives the worst correlation with 
experimental results among the three 
methods evaluated. This occurs, in 
part, because the theory does not 
account for the influence of the time- 
average loading. As discussed in the 
preceding section, it is unclear why 
this method predicts that the loading 
coefficient decreases with decreasing J, 
rather than predicting that it is 
insensitive to J. 


Near design J, the method of 
Tsakonas et al (6,7) predicts kK (py)1 
in longitudinal wakes which is close to 
the predictions by the other calculation 
methods, and reasonably close to 
experimental results. In contrast, in 
inclined flow the method of Tsakonas et 
al predicts much smaller values of 
K(Fx)1 than either the other calcula- 
tion methods or the experiments even at 
design J; see Figure 9. In particular, 
near design J the predictions of k(py)1 
by Tsakonas et al agree with the predic-—- 
tions of Kerwin and Lee in longitudinal 
flow but not in inclined flow. This sug- 
gests that the tangential inflow veloci- 
ties significantly influence the 
periodic blade loads in a manner which is 
considered by the method of Kerwin and 
Lee but which is not properly considered 
by the method of Tsakonas et al. 


*Recall that Propeller 4661 rotates left- 
hand and propeller 4710 rotates right- 
hand, and the phase angle is the angle 
of the blade reference line which the 
maximum positive loading occurs measured 
from the upward vertical positive in the 
direction of propeller rotation. 


The unsteady method of Tsakonas et a 
al predicts that the maximum value of 
(Fx)1 in the longitudinal wake will 
occur at approximately 10 degrees after ine 
the experimental value for Propeller 
4661 and approximately 10 degrees before 
the experimental value for Propeller 
4710; see Figure 13. 


A moderately-loaded unsteady 
lifting surface theory developed by 
Kerwin and Lee at MIT. 

A refinement by Kerwin to the 
method of Kerwin and Lee which 
considers the inclination of the 
propeller slipstream for opera- 
tion in inclined flow. 


The variations of kK(px)1_ The results are as follows: 
predicted by the method of Kerwin and iL, In inclined flow, all four of the 
Lee (8) agrees more closely with the calculation methods evaluated 
experimental results than either of the consistently underpredicted the 
other two methods evaluated. However, experimental values of the periodic 
the agreement with experimental results propeller blade loads. The method 
is not as good as that obtained in ~ of Kerwin, which considers the 
inclined flow by the method of Kerwin. inclination of the slipstream 


relative to the propeller axis, 
produced the best correlation with 
experimental values, and the method 
of Tsakonas et al produced the 
worst correlation with experimental 
values. The correlation of each 
method evaluated is as follows: 


The magnitude and the variation of 
the phases with advance coefficient 
predicted by the method of Kerwin and 
Lee agree with the experimental results 
substantially better than do those 
predicted by either of the other two 
predictions evaluated. Except for J> 
1.0 on Propeller 4710, the phases a. 
predicted by Kerwin and Lee are within 5 
degrees of the experimental phases. 


The method of Tsakonas et al 
predicts the first harmonic 
loading to be approximately 35 
percent of the experimental 
value at the lowest advance 
coefficient investigated, 65 
percent of the experimental 
value at the design condition 


In summary, all three of the 
calculation procedures predicted 
K(px)1 to within 20 percent of the 
experimental values at design J. 


However, the agreement at substantially 
off-design J was not as good. In 
general, the method of Kerwin and Lee 
gave the best agreement with experimen- 
tal results in the magnitude and trend 
of k(px)1 Over a range of J. In addi- 
tion, the method of Kerwin and Lee gave 
the best agreement with experimental 
results in phases and in the trends of 
the variation of phases with J. 


SUMMARY AND CONCLUSIONS 


An investigation was undertaken to 
evaluate the periodic single-blade loads 
on propellers and the influence of a 
nearby solid boundary on the periodic 
blade loads in inclined flow and in 
circumferentially non-uniform 
longitudinal velocity fields. 
Systematic model experiments were 
conducted and the results were 
correlated with predictions by the 
following methods: 


and 80 percent of the 
experimental value at the 
highest advance coefficient 
investigated. This method 


predicts a variation of first 
harmonic loading coefficients 
with advance coefficient which 


is substantially different 


from that obtained experiment- 


ally. This method predicts 
that the maximum values of 
loading will occur at 


approximately 15 to 20 degrees 


of blade angular position 
before the experimentally 
determined angle. 


b. The quasi-steady method of 
McCarthy underestimates the 


first harmonic loading by 10 
to 30 percent over a range of 
advance of coefficients with 


closer agreement at higher 
values of J. This method 


predicted a similar trend to 
the experimental data; i.e., 
that the loading coefficient 
decreases with increasing J. 
This method predicts that the 
maximum values of loading will 
occur at approximately 10 to 
20 degrees before the experi- 
mentally determined angle. 


dbo A simple quasi-steady procedure 
developed by McCarthy at DTNSRDC 
which utilizes the open water 
characteristics of the propeller. 

Bc A linearized lightly-loaded 
unsteady lifting surface theory 
developed by Tsakonas and his 
colleagues at Davidson Laboratory. 
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ce The unsteady theory of Kerwin 
and Lee showed, in general, 
somewhat better agreement with 
experimental data over a range 
of advance coefficient than 
either the quasi-steady method 
of McCarthy or the unsteady 
method of Tsakonas et al. 
This applies to both ampli- 
tudes and phases. This method 
underpredicted the periodic 
loads by approximately 5 to 20 
percent of the experimental 
values, with closer agreement 
being obtained for the Fy 
and My components at the 
higher values of J. This 
method predicts the phases to 
within 5 to 15 degrees of the 
experimental values for most, 
but not all, conditions 
evaluated. In general, the 
agreement is better for the 
higher values of advance 
coefficient J; i.e., J near or 
higher than the design value. 


ais The method of Kerwin, which is 
the only method evaluated that 
considers the inclination of 
the slipstream relative to the 
propeller axis, produced 
consistently the best agree- 
ment with experimental data, 
including both amplitudes and 
phases over the range of 
parameters evaluated. In 
particular, for low values of 
J this method yielded substan- 
tially better correlation with 
experimental results than any 
of the other methods evaluated. 
At high values of advance 
coefficient this method yields 
substantially the same results 
as the method of Kerwin and 
Lee. Thus, the inclination of 
the propeller slipstream rela- 
tive to the propeller axis can 
significantly influence the 
periodic propeller blade 
loads, and the importance of 
this inclination increases 
with increasing time-average 
loading. 


In longitudinal flow, all of the 
calculation procedures evaluated 
predict the amplitude of the 
periodic single-blade axial force 
to within 20 percent of the experi- 
mental values at design J. However, 
the agreement at substantially 
off-design J is not as good. In 
general, the method of Kerwin and 
Lee agrees best with experimental 
results in the magnitude and the 
trend of periodic single-blade 


axial force coefficient over a 
range of J. In addition, the 
method of Kerwin and Lee agrees the 
best with experimental phases and 
in the trends of the variation of 
phases with J. 


36 The periodic blade loads in 
inclined flow are not significantly 
influenced by the presence of a 
nearby boundary. 
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1 SILVERMAN 
U MARYLAND/GLEN MARTIN INST 
CIT 
1 AERO LIB 
1 ACOSTA 
1 PLESSET 
1 WU 
CATHOLIC U 
COLORADO STATE U/ALBERTSON 
U CONNECTICUT/SCOTTRON 
CORNELL U/SEARS 


FLORIDA ATLANTIC U OE LIB 


Copies 


HARVARD U 
1 MCKAY LIB 
1 BIRKOFF 
1 CARRIER 


U HAWAII/BRETSCHNEIDER 
U ILLINOIS/ROBERTSON 


U IOWA 
1 TIHR/KENNEDY 
1 IHR/LANDWEBER 


JOHNS HOPKINS U 
1 PHILLIPS 
1 INST COOP RES 


U KANSAS CIV ENGR LIB 
KANSAS ST U ENGR EXP/LIB 
LEHIGH U FRITZ ENGR LAB LIB 
LONG ISLAND U 


ORI, INC 
2 KIM 
1 SCHNEIDER 
1 WILLIAMS 


MIT 
1 BARKER ENGR LIB 
2 OCEAN ENGR/KERWIN 
1 OCEAN ENGR/LEEHEY 
1 OCEAN ENGR/NEWMAN 


U MINNESOTA SAFHL 
1 KILLEN ; 
1 SONG 
1 WETZEL 


STATE U MARITIME COLL 
S U ARL LIB 
1 ENGR DEPT 
1 INST MATH SCI 


NOTRE DAME ENGR LIB 


Copies 


PENN STATE U ARL 
LIB 
HENDERSON 
GEARHART 
PARKIN 
THOMPSON 


| el 


PRR Ee 


PRINCETON U/MELLOR 
RENSSELAER/DEPT MATH 
ST JOHNS U 
SWRI 

1 APPLIED MECH REVIEW 

1 ABRAMSON 

1 BURNSIDE 
BOEING ADV AMR SYS DIV 
BOLT BERANEK AND NEWMAN 

1 BROWN 

1 JACKSON 
BREWER ENGR LAB 
CAMBRIDGE ACOUS/JUNGER 
CALSPAN, INC/RITTER 
STANDFORD U/ ASHLEY 
STANDFORD RES INST LIB 
SIT DAVIDSON LAB 

1 erenB 

1 BRESLIN 

1 TSAKONAS 
TEXAS U ARL LIB 
UTAH STATE U/JEPPSON 
WEBB INST 

1 WARD 

1 HADLER 


WHOI OCEAN ENGR DEPT 
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Copies 


WPI ALDEN HYDR LAB LIB 
ASME/RES COMM INFO 
ASNE 
SNAME 
AERO JET-GENERAL/BECKWITH 
ALLIS CHALMERS, YORK, PA 
AVCO LYCOMING 
BAKER MANUFACTURING 
BATH IRON WORKS CORP 

1 HANSEN 

1 FFG PROJECT OFFICE 
BETHLEHEM STEEL SPARROWS 
BIRD-JOHNSON CO 

1 CASE 

1 RIDLEY 

1 NORTON 
DOUGLAS AIRCRAFT 

1 TECHNICAL LIBRARY 

1 SMITH 
EXXON RES DIV 

iL Lat 

1 FITZGERALD 
FRIEDE & GOLDMAN/MICHEL 


GEN DYN CONVAIR 
ASW-MARINE SCIENCES 


GIBBS & COX 
1 TECH LIB 
1 OLSON 
1 CAPT NELSON 


GRUMMAN AEROSPACE/ CARL 


HYDRONAUTICS 
1 ETTER 
1 SCHERRER 
1 LIBRARY 


Copies Copies 
iL INGALLS SHIPBUILDING 1 TRACOR 
i INST FOR DEFENSE ANAL il UA HAMILTON STANDARD/ CORNELL 
il ITEK VIDYA 


CENTER DISTRIBUTION 
a LIPS, INC 


Copies Code Name 
1) OLE RreETON) Rie EN. 
TLETON GR CORP/REED i te en 
1 LITTON INDUSTRTES 1 1102.1 Nakonechny 
1 LOCKHEED/WAID 1 15 Morgan 
1 MARITECH, INC/VASSILOPOULIS 7 1202 Poe ee 
1 152 Lin 
3. HYDRODYNAMICS. RESEARCH ASSOCIATES, : Al een 
INC 
COX iL 1528 Lee 
1 VALENTINE i 1528 Reed 
1 NELKA 1 1528 Dobay 
aL 154 McCarthy 
1 NATIONAL STEEL & SHIPBLD 
a e * EE be 1 1543 Cumming 
1 NEWPORT NEWS SHIPBLDG LIB iL 1544 Brockett 
30 1544 Boswell 
1 NIELSEN ENGR/SPANGLER 10 1544 Jessup 
1 HYDROMECHANICS, INC/KAPLIN . a Huang 
1 1556 Santore 
1 NAR SPACE/UJIHARA iL 1556 Wisler 
alt 1556 Jeffers 
1 PROPULSION DYNAMICS, INC 
il, L/P Krenzke 
1 PROPULSION SYSTEMS, INC i pl bie Secale 


1 SCIENCE APPLICATIONS, INC/STERN 


1 19 Sevick 
1 GEORGE G. SHARP 1 19 Strasberg 
1 1903 Chertock 
1 SPERRY SYS MGMT LIB/SHAPIRO 1 lace Oey 
2 SUN SHIPBLDG 1 1962 Noonan 
1 LIBRARY 
1 NEILSON 1 2814 Czyryca 
1 ROBERT TAGGART 10 5211.1 Reports Distribution 
iL SP Al Library (C) 
iL TETRA TECH PASADENA/CHAPKIS 1 BOD 2 ‘ieiarecren (OS 
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DTNSRDC ISSUES THREE TYPES OF REPORTS 


1. DTNSRDC REPORTS, A FORMAL SERIES, CONTAIN INFORMATION OF PERMANENT TECH- 
NICAL VALUE. THEY CARRY A CONSECUTIVE NUMERICAL IDENTIFICATION REGARDLESS OF 
THEIR CLASSIFICATION OR THE ORIGINATING DEPARTMENT. 


2. DEPARTMENTAL REPORTS, A SEMIFORMAL SERIES, CONTAIN INFORMATION OF A PRELIM- 
INARY, TEMPORARY, OR PROPRIETARY NATURE OR OF LIMITED INTEREST OR SIGNIFICANCE. 
THEY CARRY A DEPARTMENTAL ALPHANUMERICAL IDENTIFICATION. 


3. TECHNICAL MEMORANDA, AN INFORMAL SERIES, CONTAIN TECHNICAL DOCUMENTATION 
OF LIMITED USE AND INTEREST. THEY ARE PRIMARILY WORKING PAPERS INTENDED FOR IN- 
TERNAL USE. THEY CARRY AN IDENTIFYING NUMBER WHICH INDICATES THEIR TYPE AND THE 
NUMERICAL CODE OF THE ORIGINATING DEPARTMENT. ANY DISTRIBUTION OUTSIDE DTNSRDC 
MUST BE APPROVED BY THE HEAD OF THE ORIGINATING DEPARTMENT ON A CASE-BY-CASE 
BASIS. 
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